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The critical behavior of the 1/5-depleted square- lattice Ising model with nearest neighbor fer- 
romagnetic interaction has been investigated by means of both an exact solution and a high- 
temperature series expansion study of the zero-field susceptibility. For the exact solution we employ 
a decoration transformation followed by a mapping to a staggered 8-vertex model. This yields a 
quartic equation for the critical coupling giving K c (= fiJ c ) = 0.695. The series expansion for the 
susceptibility, to 0(K 18 ), when analyzed via standard Pade approximant methods gives an esti- 
mate of K c , consistent with the exact solution result to at least four significant figures. The series 
expansion is also analyzed for the leading amplitude and subdominant terms. 
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I. INTRODUCTION 

Exact solutions of lattice models play an important 
role in the study of phase transition and critical phenom- 
ena. In his seminal work, Onsager solved the two dimen- 
sional (2D) square lattice Ising model (2D- Ising) exactly 
PP. Solutions have been obtained for other regular 2D 
lattices [2J. A number of complex configurations such 
as the Union Jack, the bathroom tile (or 4-8), the 4-6, 
and the l/9 th depleted lattice models have also been in- 
vestigated [3HH]- The 1/5-th depleted antiferromagnetic 
S=l/2 Hcisenberg model has been an earlier topic of in- 
vestigation [9] . In this article we predict the critical point 
of a previously unexplored 1/5-th depleted Ising model 
on a square lattice. The Ising model Hamiltonian is given 
by 

u = -jy~i - hy^sj, (i) 

(hi) ' 

where Sj is the classical dimensionless Ising variable at 
site i taking the values ±1. J > (ferromagnetic interac- 
tion) and h (magnetic field) denote constant parameters 
with dimensions of energy. The structure of the depleted 
lattice with every 5 th missing site is shown in Fig. Ill 

In principle all planar Ising models (i.e. with non- 
crossing bonds) are solvable by the Pfaffian method [2J. 
The method has been utilized to solve a variety of lattice 
models including the 8-vertex model. The 8-vertex model 
has been investigated both for translationally invariant 
and staggered vertex weight [TO] E]. In the staggered 
model the vertex weights are allowed to vary taking dif- 
ferent values on the staggered plaquettes of the square 
lattice. The relevance of the staggered model lies in its 
relationship to a number of important models in statis- 
tical mechanics - the percolation model [TJ], the Potts 
model [13], and the Ashkin- Teller model [Ml ITS] . 
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FIG. 1. (Color online) 1/5-th depleted lattice Ising model. 
Filled red circles denote points on the original lattice. Yel- 
low squares denote points on the decorated lattice. Solid blue 
lines denote bonds on the original lattice. Dashed black lines 
denote bonds on the decorated lattice. A and B represent the 
two types of plaquettes created after the decoration transfor- 
mation. The transformed lattice has a staggered configuration 
of A and B plaquettes. 



The partition sum of the 2D-Ising model has been eval- 
uated using a variety of techniques [REF]. We compute 
the partition sum to obtain the critical point of the 1 /5-th 
depleted Ising model with two methods. First, we carry 
out an exact solution by using a decoration transforma- 
tion [16] followed by a mapping to a staggered 8-vertex 
model [HlUTT] . Second, we obtain the high temperature 
series expansion (HTSE) for the zero-field susceptibility 
upto 0(K 18 ) where K = pj (j3 = l/k B T). T is the tem- 
perature and ks the Boltzmann constant. Using Pade 
approximants (PA) we analyze the series for its leading 
amplitude and subdominant terms. 

This paper is organized as follows. In Section [T] we in- 
troduce the 1/5-th depleted lattice Ising model. In Sec- 
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FIG. 2. (Color online) Ising spins on the original lattice 
are denoted by Si, i=l,2,3,4 (filled red circles). Spins on the 
transformed decorated lattice are given by cr;, i=l,2,3,4 (yel- 
low squares). Interactions between s-spins is given by K. In- 
teractions between s and cr-spins is given by K . 



tion |n] we obtain the exact solution. In Section Mil we 
compute the HTSE for the zero-field susceptibility and 
display the series coefficient results. In Section[IV]we an- 
alyze the leading and subleading amplitudes of our series 
expansion result. Finally in Section[V]we summarize and 
conclude the main results of the paper. 



II. EXACT SOLUTION 

We carry out a two step procedure to obtain the exact 
solution. In the first step a decoration transformation is 
performed. A new spin a (filled yellow squares in Fig. [I]) 
is introduced in the original lattice placed at the mid- 
point of bonds which do not lie in the small square (see 
Fig.[lJ [16 . The interactions on these bonds are replaced 
by new interactions with the a spins to obtain the rela- 



tion between the spins (s) in the original lattice and the 
decorated lattice (cr). The scaled interactions K and K 
are defined by 



(2) 



Carrying out the summation we obtain the following re- 
lation between K and K 



„2A" 



cosh(2A' ). 



(3) 



The decoration process introduces two different types of 
plaquettes A and B as shown in Fig. [2] We map this 
configuration to a staggered 8-vertex model. This trans- 
forms the 1 /5-th depleted Ising model to a square lattice 
with more complex interactions with nearest neighbour 
coupling K, a diagonal next-nearest-neighbour coupling 
L, and a four-spin coupling M [17] . The transformed Z 
is 

Z(oi, cr 2 ,cr 3 , C4) = ^exp[K (si<7i + S 2 cr 2 + .S 3 cr 3 + S4CT4)] 

W 

x exp[A'(sis 2 + s 2 s 3 + S3S4 + S1S1)]. 

, . (4) 

Enumerating the summation over the original lattice 
spins, Si (i=l,2,3,4), we can write 

Z(cr 1 ,(T2, 0-3,0-4) = A(K) exp[A v (cria 2 + 020-3 + <t 3 <t 4 + 01O4)] 
x exp[L(oio 3 + 0-20-4)] exp[M(o-io- 2 o-30- 4 )], 

(5) 



with 

Ae 4K+2L+M = Ae -2K' + ^2K + ^K-IK + & 4k' ' +4K + 2& -4K + At = p^K), 
Ae -4K+2L+M = 4e -2A-' + 4e 2A-' + & -i K ' -iK + £ 4 K ' -4K + 2& 4K + 4 = P ^ K ) } 
Ae -2L+A/ = e -4K + Ae -2K + ^2K + e 4K + 2& -4K + 2& 4K + 3 = p^K), 
Ae -M = e -4K + 2e -2K' + 2e 2K + + e -2K -4K + £ 4K-2k' + & 2K ' -AX + £ 2K +4K + g = 



The vertex weights satisfy the free fermion condition |llj . 
Using Eqs. 13 and 21 from Ref. [IT] the condition for the 
critical point for our model is 

Pi{K) + P 2 (K) - 2P 3 {K) - m(K) = 0. (10) 

The above equation can be written in terms of the vari- 
able x — e 2K as 

x 4 - 4x 3 -1 = 0. (11) 

Solving the above equation we obtain a physical root at 
x — 4.015445, yielding a critical coupling value of K c — 



(6) 
(7) 
(8) 

P 4 (K). (9) 
I 

0.695074. In terms of the variable v = tanh(A') Eq. [II 
takes the form 

v 4 + 4v 3 -l = 0, (12) 

giving v c =0. 601232 as the solution of the critical point. 
III. HIGH TEMPERATURE SERIES 

The HTSE technique is one of the most effective ap- 
proaches to study critical phenomena |18j . Much work 
has been devoted to the HTSE of the Ising model [T51 - 
I21j . Thermodynamic properties are derivable from the 
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TABLE I. High-temperature series expansion coefficients, 
a r , of the zero-fieid susceptibility for the 1/5 depleted fsing 
model. The expansion parameter is v = tanh(/9J). 



Order 


Coefficient 


1 


0.300000000000D+01 


2 


0.600000000000D+01 


3 


0.120000000000D+02 


4 


0.220000000000D+02 


5 


0.400000000000D+02 


6 


0.740000000000D+02 


7 


0.136000000000D+03 


8 


0.246000000000D+03 


9 


0.444000000000D+03 


10 


0.782000000000D+03 


11 


0.137200000000D+04 


12 


0.240600000000D+04 


13 


0.420800000000D+04 


14 


0.738600000000D+04 


15 


0.129240000000D+05 




U.zzoy4UUUUUUUiJ-t-U0 


17 


0.387280000000D+05 


18 


0.667820000000D+05 


partition function 




Z = ^cxp(-/?H), 

S 


(13) 




4-/3^5,), (14) 

i 


The zero field susceptibility, is . 


given by 



Using the identity below for both the regular and the 
field term 

exp(Ksi-Sj) — cosh if (1 + vsiSj), (16) 

we can construct a graphical expansion. Each bond car- 
ries a factor of ws^Sj and, in addition, each site has a 
factor of either 1 or TSfe. Only those graphs with pre- 
cisely two factors of TSk contribute to the above equa- 
tions. As a result the graphs which contribute are those 
with precisely two vertices of odd degree, those to be 
compenstated by the two TSk factors. We then obtain 
the following result 

(coshlfF^coshW = § EK 1 + ™<*i) IB 1 + "*) 

(17) 

where v = tanh j3 J and r = tanh f3h. The high tempera- 
ture susceptibility can be expanded in the form 

oo 

p- l X {v) = l + Y,arV r (18) 

r=l 



The coefficients be related to the graph counting 

problem and evaluated exactly [18] . The computed series 
expansion coefficients for the zero-field susceptibility of 
the 1/5-th depleted Ising model are listed in Table [i] 
IV. SUSCEPTIBILITY ANALYSIS 

The universality hypothesis in critical phenomena im- 
plies that thermodynamic quantities are not sensitive to 
the microscopic details of a system near a critical point 
[22l [23] . It is known from earlier work that near the tran- 
sition point the high temperature susceptibility, xi v ) i of 
the 2D-Ising model on all 2D lattices has an asymptotic 
form. For our model we can express the susceptibility as 

*v)=A.(l--) +-(19) 

with v c —0. 601232. To analyze the x( v ) series for its pole 
and its leading and subleading amplitude we first con- 
sider constructing the series 

fi(v) = lx(v)} 4/7 ~ A 4 J 7 (1 - v/vc)- 1 + ■■■ (20) 

Direct PA's to fi(v) give a consistent pole at v c ~ 0.6015 
± 0.0002. This result is close to the exact solution value 
of v c =0. 601232. The residues, which are estimates of 
P^J 1 v c are all in the range 0.486 - 0.490. Considering 
the value to be 0.488 we obtain A Q ~ 0.694. A more 
consistent set of results can be obtained by constructing 
the series 

f2{v) = (l-v/v c f/\{ V ) 

^A Q + terms which vanish at v c (21) 

and forming PA's to /2(f)- Evaluating these at 
v c =0.601232 gives A ~ 0.687 ± 0.001. To obtain the 
subdominant contribution we analyze the function 

h(v) = (l-v/v c f/* X (v)f 2 (v) 
~Ai + terms which vanish at v c . (22) 

PA's of fs(v) computed at v c —0. 601232 provide a consis- 
tent set of estimates for Ai ~ 0.708 ± 0.001. 

The above analysis can be repeated with the HTSE \ 
series expressed in the K variable, 

/ K \- 7/i ( K\- 3/4 

(23) 



TABLE II. Pade approximation analysis of the high- 
temperature series expansion coefficients of the zero field sus- 
ceptibility for the 1/5 depleted Ising model. The critical cou- 
pling constant and the leading amplitude is listed below. 



(N,D) 


Vc 


Ac 


Ai 


(10,8) 


0.601461 


0.686077 


0.707575 


(9,9) 


0.601405 


0.686598 


0.707919 


(8,10) 


0.601437 


0.686265 


0.707625 


(9,8) 


0.601554 


0.686771 


0.707614 


(8,9) 


0.600934 


0.686810 


0.707385 
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To obtain a consistent set of critical coupling value, K c , 
we perform the PA analysis on a x 4 ^ 7 series. The leading 
amplitude, C D , can be computed by investigating the PA 
analysis of (1 — K/K c ) 7 ^ 4 \- However, such an analysis 
does not lead to a consistent set of values for C\. We 
therefore obtain both Co and C\ from Aq and A\, To 
do so we expand the x( v ) series in a Taylor series in 1- 
K/K c up to 2nd order to obtain the following 

K c = 0.695 ±0.001, (24) 
C* = 1.167 ± 0.001, Ci = 0.036 ± 0.001. (25) 

The agreement of the series estimate of K c with our result 
from Eq . |11| provides confirmation that our analysis in 
Section Ml is correct. 

V. CONCLUSION 

We have obtained the critical point exactly and the 
estimated values of the leading two amplitudes of the 
asymptotic form of the zero-field susceptibility for the 
Ising model on an unusual lattice obtained by regularly 
removing the 1/5-th of the sites of a square lattice. There 
has not been to our knowledge, any systematic study of 
Ising models on depleted lattices. It is worth noting, 
however, that the familiar honeycomb and kagome lat- 
tices result from particular 1/3-rd and 1/4-th depletions 
of the triangular lattice. The 1/5-th depleted square lat- 
tice, considered here, is in fact realized in the material 
CaV^g, but not as an Ising system. 

To obtain the critical point, we relate partition func- 
tion of the model to that of a staggered eight-vertex 
model and use established results for that model. The 
vertex weights satisfy a 'free fermion' condition, confirm- 
ing that the model lies in the normal Ising universality 
class. We note that in the 1/5-th depleted lattice there 
are two classes of nearest-neighbor bonds. While we have 
only considered the case of equal strengths, our transfor- 
mation method applies equally well to the more general 



case of different couplings J, J . 

We have also derived an 18-term high temperature se- 
ries for the zero-field susceptibility. Because of the low co- 
ordination number and open structure of the lattice, the 
number of graphs that contribute at a given order is much 
reduced. As a consequence, the series is not as well be- 
haved as that of the parent square lattice. However, using 
standard PA methods, we obtain an estimate of the criti- 
cal temperature in good agreement with the exact value. 



As explained in Section IV we obtain rather precise es- 



timate of the leading two amplitudes in the asymptotic 
form of x( v ) near the critical point and, rather less pre- 
cise, estimates of the amplitudes in the K-representation. 

For the square lattice exact expressions of the spin-spin 
correlation functions allow these amplitudes to be ob- 
tained, essentially exactly, from the solution of a Painlcve 
equation III 24J. In addition an exact result relat- 
ing the coefficients Co, Ci, viz., Ci/Cq =V^K c /8, has 
been proven. It is not clear whether a similar calculation 
could be done for the depleted lattice. However, the ratio 
Ci/Co, in this case, does not appear to satisfy a simple 
relationship of the above type. 

Finally, we remark that depletion of any regular lat- 
tice will reduce the average coordination number and this 
leads to a less rigid structure. Hence the ordered state 
will be less robust to thermal fluctuations, and the criti- 
cal temperature will be lowered. This is seen in our study, 
with kBT c /J being reduced by some 36%, from 2.2692. . . 
to 1.4387. . . . 
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